Draft version April 27, 2010 

Preprint typeset using WT^i. style cmulatcapj v. 08/13/06 



O 

5-1 

< 



AN ANALYTIC MODEL OF THE PHYSICAL PROPERTIES OF GALAXY CLUSTERS 

G. Esra Bulbul 1 , Nicole Hasler 1 , Massimiliano Bonamente 1 ' 2 and Marshall Joy 2 

Draft version April 27, 2010 

ABSTRACT 

We introduce an analytic model of the diffuse intergalactic medium in galaxy clusters based on 
a polytropic equation of state for the gas in hydrostatic equilibrium with the cluster gravitational 
potential. This model is directly applicable to the analysis of X-ray and Sunyaev-Zeldovich Effect 
observations from the cluster core to the virial radius, with 5 global parameters and 3 parameters 
describing the cluster core. We validate the model using Chandra X-ray observations of two poly- 
tropic clusters, MS 1137.5+6625 and CL J1226. 9+3332, and two cool core clusters, Abell 1835 and 
Abell 2204. We show that the model accurately describes the spatially resolved spectroscopic and 
imaging data, including the cluster core region where significant cooling of the plasma is observed. 
Subject headings: X-rays: galaxies: clusters-galaxies: individual (MS 1137.5+6625, CL J1226. 9+3332, 
Abell 1835, Abell 2204) 
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1. INTRODUCTION 

Galaxy cluster masses play an important role in ad- 
dressing fundamental physical and cosmological prob- 
lems, such as the measurement of the clus t er ga s 
mass fraction (jAllen et all 120081: lEttori et all 12 009). 
the e v olution of the growth of str ucture (|Mantz et al.l 
120081 120091 : IVikhlinin etall 12009ft. and the gravita- 
tional sedimentation of ions (IChuzhov fc Nusserl 120031 : 
iPeng fc Nagal [200l iShtvkovskiv fc Gilfanovl ^OlOf TA 
vital tool for the measurement of cluster masses is 
the diffuse hot intergalactic medium which can be 
detected prim arily through its bright X-ray emission 
(jSarazinl I1988T) . or through the Sunyaev-Zeldovich Ef- 
fect (SZE, ISunvaev fc Zel'dovichl [1973 ICarlstrom etall 
l2002h . A variety of models are used to de- 
scribe the distributio n of the gas, from the simple 
isothermal model (|Cavaliere fc Fusco-Femianol 119761 : 
Birkinsha w et al.l 119911 ) to more complex models that 
describe eith er the X-ray properties (gas density and 
temperature, IVikhlinin et ; al.ll200d ICavaliere et alJl2009ft 
or the SZE properties (sa,s P ressure - iNagai et al.l [20071 : 



iMroczkowski et alj|2009t Urnaud et afl feoOQ) 



We investigate a model of galaxy clusters based on 
an analytic distribution for the c luster mass density in- 
spired by the lNavarro et al.l (119961) dist ri bution, which we 
genera lize following iSuto et al.l (119981 . lAscasibar et al.l 
(2003|) and lAscasibar fc Diegol (2008) to include a vari- 
able asymptotic slope at large radii. This mass density is 
combined with a polytropic equation of state for the gas, 
to provide self-consistent density, temperature and pres- 
sure profiles for a plasma in hydrostatic equilibrium. The 
use of a polytropic equation of state for the cluster gas 
has al so be en recently propo sed by lAscasibar fc Diegol 
(2008) and iBode et all (120091) and tested observation- 



ally by ISanderson fc Ponmanl (|2009l ) . 

In this paper we derive analytic radial profiles for 
the physical quantities (temperature, density and pres- 
sure) relevant to X-ray and SZE observations, and 
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present an application of these models to high res- 
olution Chandra observations of the galaxy clusters 
MS 1137.5+6625, CL J1226.9+3332, Abell 1835 and 
Abell 2204. Applications of this new model include 
measurement of gas mass fraction from joint X-ray 
and S unyaev-Zeldovich Effect observations ([Hasler et al.l 
|2010( ) and the effect of He sedim entation on X-ray 
mass estimates (Bul bul et aJ]l2010h . This paper is or- 
ganized as follows: in [21 we describe our model, in fj3] 
we present the application of the model to Chandra X- 
ray observations of MS 1137.5+6625, CL J1226.9+3332, 
Abell 2204 and Abell 1835, and in J|]we perform a com- 
p arison between our mass m easurements and the results 
of Mroczkowski et al. (2009). In fj5]we present our con- 
clusions. In the analysis of the Chandra data we assume 
the cosmological parameters h = 0.73, £Im = 0.27 and 
n A = 0.73. 

2. A MODEL OF THE INTERGALACTIC MEDIUM BASED 
ON HYDROSTATIC EQUILIBRIUM AND THE 
POLYTROPIC EQUATION OF STATE 

2.1. The Mass Density Distribution 

The cluster gravitational potential is dominated 
by dark matter, with the intergalactic medium and 
stars contr ibuting less than approximately ~20% of 
the mass (I Allen et all I2004L I2008L lEttori et al.l 120091 : 
Vikh linin et al.ll2009D . We therefore start with a total 
mass density distribution that is o btained as a general- 
ization of the [Navixro^Eid. (1996) profile: 



Ptot{r) 



(r/r s )(l + r/r s )P 



(1) 



where pi is the normalization constant, r s is a charac- 
teristic scale radius and (3+1 is the slope of the density 
distribution at large radii. Equation [1] is a si mplified ver- 
sion o f the density distribution introduced bv lSuto et al.l 
(fl998h . 

The total mass enclosed within radius r can be found 
by taking the volume integral of the density (Equation[lJ : 



M(r) 



An pi 



1/(1 - P) - r/r s 



09-2) Vj" - 1 (1 + r/r,)"- 



(2) 
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Equation [2] is indeterminate at j3 = 2; the limit at 
/3 = 2 can be determined using L'Hospital's rule: 



M(r) = 47rp,rf ln(l + r/r s ) 



r/r s 



1 + r/r s 



\J3 = 2] (3) 



and therefore the mass is a continuous function of /3 
with no discontinuity at /3 = 2. 

The gravitational potential at a distance r is found by 



d(j){r) = GM(r)/r 2 dr 



(4) 



using the boundary condition 4>(oo) — 0. Equation U 
can be integrated analytically 



(r) = 0o 



where 



1 (l + r/r s )' 3 - 2 - 1 
(0-2) r/r s (l + r/r s )^ 2 



(5) 



(6) 



The potential at /3 = 2 is also found by using 
L'Hospital's rule, 

H r)-*(*i±&l) (7, 

Therefore the gravitational potential is a continuous 
function of /3 with no discontinuity at /3 = 2. Figure 
Q] shows the radial distribution of the gravitational po- 
tential for 1.0 < f3 < 3.0. The limiting value of j3 = 1 is 




Fig. 1. — Normalized gravitational potential for various values 
of the parameter. 



shown in Figure [1] corresponding to a constant potential. 

2.2. Gas Density and Temperature Profile 

The diffuse gas is assumed to be in hydrostatic equilib- 
rium with the gravitational potential. Assuming spheri- 
cal symmetry, 



1 



dP P 



pm p n e (r) dr 



dcj){r) 
dr 



(8) 



where P e (r) — n e (r)kT(r) is the electron pressure, 
G denotes the gravitational constant, m p is the proton 
mass, p is the mean molecular weight of the plasma, k is 



the Boltzmann constant and n e (r) is the electron num- 
ber density. In order to solve Equation[8j we assume that 
the gas follows a polytropic equation of state, 



,poly( r ) 



Tpoly (^*) 

T 



(9) 



where n is the polytropic index, n e Q and To are the 
values of the number density and temperature at r = 0. 
The polytropic index n > is a free parameter of the 
model, with the li mit n — » oo descr ibing an isothermal 
distribution of gas (Eddingto~nl ll926(h 

The temperature profile is obtained as a function of 
the gravitational potential from Equations [8] and El 



Tpoly(r^ - 



firrip 



1) k 

and therefore, using Equation O 



4>(r) 



(10) 



Tpoiy(r) = T 



1 (1 + r/rs)?- 2 - 
(0-2) r/r s (l + r/r s )P- 2 J ' 



(11) 



where the normalization constant To is obtained from 
Equations [5] and [TU1 



To = 



AirGfirrip r 2 pi 
k{n+l) 08-1)' 



(12) 



Equation [TOl shows that T po i y (r) cx <fi(r), and therefore 
Figure Q] also describes T po i y (r) as function of radius. 
Equation [T2] links the gas temperature to the normaliza- 
tion of the matter density pi, and therefore the depth of 
the gravitational potential can be determined from the 
observed temperature profile. 

Using the relation between temperature and gas den- 
sity provided by the polytropic relation (Equation [§]) , 
the polytropic gas density profile is 



, 1 (l + r/r,p-l V 
n e ,poi y {r) = n e0 [ — n , „i„ \/3-2 } • ( 13 ) 



08-2) r/r s (l + r/r s y 



2.3. The Electron Gas Pressure 

The gas pressure is obtained using the ideal gas law 
T e (r) — n e ^ P oiykT po iy(r) , 



Pe{r) = PeO 



1 (1 +r/r s )^~ 2 - 1 
08-2) r/r s (l + r/r s )P- 2 



(14) 



In the limit (3 2 the pressure is analytically described 



by 



Pe{r)=PeO 



/ ln(l+r/r 8 ) 
V r/r„ 



71+1 



\P = 2\. (15) 



The electron pressure for this model has only 4 
free parameters, and it is suitable for the analysis of 
Sunyaev- Zeldovich Effect observations of galaxy clusters 
(|Hasler et al.ll2010h . 
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2.4. Cool Core Clusters 

Although the temperature profile predicted by the 
polytropic model provides a good description at inter- 
mediate to large radii, cool core clusters feature a signif- 
icant temperature drop in the central region which can- 
not be approxim ated by a polytropic equation of state 
(see for example, [y ikhlinin et al. 2005; Sanderso n et al.l 
120061: IVikhlinin etall l2006t iBaldi et al.l 120071) . For cool 
core clusters we introduce a modified temperature profile 



T(r) =T po i y (r)T coo i(r) 



(16) 



where T po i y (r) is the temperature profile according 
to the polytropic equation of state (Equation [TTj) and 
Tco oi(r) is a phenqmenolo gical core taper function used 
bv IVikhlinin et"aH (l200l : 



7~cool 



(r) 



a + (r/r coo i)' y 
1 + (r/r coo ip : 



(17) 



where < a < 1 is a free parameter that measures the 
amount of central cooling and r coo i is a characteristic 
cooling radius. The temperature profile modified by the 
core taper function is shown in Figure[2]for representative 
values of parameters r coo ;, 7 and a. 

Therefore, the explicit temperature profile for cool core 
clusters is given by 



T{t) = To 



1 



(l + r/r s ) 



0-2 



1 



T C ooi{r). (18) 



.08-2) r/r s {\ + r/r s f- 

In order to calculate the density distribution for cool 
core clusters, we assume that the pressure distribution is 
the same as in the polytropic case (Equation [TJ}. There- 
fore, the electron density is given by 



n e (r) = 



Pe(T) 

kT{r) 



1 {l + r/r s Y 3 - 2 - 1 
09-2) r/r s (l + r/r s )?- 2 



cool 



(r) 



(19) 

The behavior of the gas density for various core taper 
parameters is shown in Figure [5J 

For hydrostatic equilibrium to be satisfied, these modi- 
fied density and temperature distributions require a mod- 
ified total mass distribution: 



M(r) = 



1/(1 - /3) - r/r s 



i(r). 



(20) 



09-2) V/3 - 1 (1 + r/r^" 1 

The only difference between the cool core total mass 
distribution (Equation [2D]) and the polytropic total mass 
distribution (Equation [2]) is the term t coo i (r) , which is 
significant only at small radii. At large radii , the effect 
of the core taper vanishes, and the thermodynamics of 
the gas is described by the polytropic equation of state. 

3. APPLICATION TO Chandra OBSERVATIONS OF 
CLUSTERS 

3.1. Chandra Data Reduction and Analysis 

We use deep Chandra ACIS-I observations of four 
galaxy clusters to validate our models: two clusters which 
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Fig. 2. — Solid red lines are the normalized temperature and 
density profiles for a taper function with parameters a = 0.3, 
r cool — t b and 7 = 2.0, and variable values for 0; blue dashed 
lines are the models without the core taper function. 

do not have a cool core component, MS 1137.5+6625 and 
CL J1226.9+3332, and two cool core clusters, Abell 2204 
and Abell 1835. The observations are summarized in 
Table [TJ As part of the data reduction procedure, we 
applied afterglow, bad pixel and charge transfer ineffi- 
ciency corrections to the Level 1 event files using CIAO 
4.1 and CALDB 4.1.1. Flares in the background due to 
solar activit y are eliminated us i ng lig ht curve filtering as 
described in iMarkevitch et ail ([2003). Filtered exposure 
times are also given in Table Q] 

For the purpose of background subtraction we use 
blank-sky observations. Given that the background is 
obtained from regions of the sky that may have different 
soft X-ray fluxes than at the cluster position, we use a 
peripheral region of the ACIS-I detector to model the dif- 
ference between the blank-sky and the cluster soft fluxes. 
This step in the analysis is particularly important for 
Abell 2204, which lies in a region of significantly higher 
soft X-ray emission than the average blank-sky region. 
Spectra and images used in this paper are extracted in 
the energy band 0.7-7.0 keV, chosen to minimize the ef- 
fect of calibration uncertainties at the lowest energies, 
and the effect of the detector background at high energy. 

Spectra are extracted in concentric annuli surround- 
ing the centroid of X-ray emission after all point sources 
were removed. An optically thin plasma emission model 
(APEC in XSPEC) is used, with temperature, abun- 
dance and normalization as free parameters. The red- 
shift and Galactic Nh of the four clusters are shown in 
Table |U 
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TABLE 1 

Cluster Sample 



Cluster z Nh a Obs. ID Exposure Time 



(cm 2 ) (ksec) 

Abcll 2204 0.152 6 5.67xl0 20 7940 72.9 

Abell 1835 0.252 6 2.04xl0 20 6880 110.0 

MS 1137.5+6625 0.784 c 9.54xl0 19 536 115.5 

CL J1226.9+3332 0.888 d 1.38xl0 20 5014 32.7 

3180 31.5 



(a) Leiden/ Argentine/B onn (LAB) Surve y, see Kal berla et"aH (2005) 
(6) Str uble Rood! fl99Sf) 
(c) IDonahue et al l 11999) 
(cO IEbeling et al.l 120011 ) 



TABLE 2 

Sources of systematic error in the Chandra data 



Source of uncertainty 


Observable affected 


Fractional Error 


Background level 


Background count rate 


5% 


Spatial variations of A a ff a 


Photon count rates 


1% 


Energy calibration of A e ff b 


Temperature measurement 


5% 



(a) Reference: http://cxc.harvard.edu/cal/ 
(6) Reference: http:/ /cxc. harvard. edu/ciao4.1/why/caldb4. l.l_hrma.html 



3.2. Systematic Uncertainties in the Chandra Data 
Analysis 

We consider possible sources of systematic uncertainty 
in the Chandra data. The blank-sky background used 
in our analysis is normalized to the high-energy back- 
ground level of each cluster observation, determined from 
peripheral regions of the ACIS detector that are free of 
cluster emission (following iMarkevitch etafll2003t l. The 
primary source of uncertainty in the background subtrac- 
tion is the choice of a peripheral region as representative 
of the background at the cluster location. Due to the 
scatter in the count rate of various peripheral regions in 
each cluster observation, we estimate a ^5% uncertainty 
in the determination of the background level from these 
Chandra observations. We use this uncertainty in the 
spectral and imaging data analysis. 

Calibration of the ACIS effective area is another sig- 
nificant source of systematic uncertainty in our analysis. 
For the spectral data used for measuring the gas temper- 
ature, the primary source of uncertainty is the low-energy 
calibration of the effective area and the presence of a con- 
taminant on the optical filter of the ACIS detector. We 
use the Chandra calibration available in CALDB 4.1.1, 
which includes a significant change in the effective area 
calibration which improves the agreement between clus- 
ters temperatures obtained with ACIS-I, and also with 
other instruments (such as XMM-Newton's EPIC). With 
this calibration of the Chandra efficiency we estimate 
that any residual systematic error in the measurement of 
cluster temperatures is of order ~5%, and add this error 
to the temperature measured in each bin. 

For the imaging data, spatially-dependent non- 
uniformities in the ACIS efficiency are a relevant source 
of possible systematic error because of the extended na- 
ture of the sources. The absolute calibration of the ACIS 
efficiency is currently at the level of 3%, with possible 
spatial variations on arcmin scales at the level of ~1%; 



we use a 1% error as additional uncertainty in the count 
rates for each annulus. 

Table [5] provides a summary of the uncertainties in- 
cluded in our analysis of the Chandra data, and refer- 
ences to the Chandra calibration information. 

3.3. Result of Model Fits and Mass Measurements 

The radial profiles of the X-ray surface brightness and 
temperature observed from the Chandra data are used 
to determine the best-fit parameters and the goodness of 
fit for MS 1137.5+6625, CL J1226.9+3332, Abell 1835 
and Abell 2204. The X-ray surface brightness is 



1 



n 2 e A ee {T)dl, 



47r(l + z) 3 

where S x is in detector units (counts cm -2 



(21) 



arcmin 

s _1 ), z is the cluster redshift, A ee (T) is the plasma emis- 
sivity in detector units (cou nts cm 3 s" 1 ) which we cal- 
culate using the APEC code ([Smith et al.ll2001l) and I is 
the distance along the line of sight. 

We validate the model using two polytropic clus- 
ters which do not have a cool core component, 
MS 1137.5+6625 and CL J1226.9+3332, and two cool 
core clusters Abell 1835 and Abell 2204. We use 
the Monte Carlo Ma rkov Chain code described in 
iBonamente et al. (2004) for the fit. The model described 
in § [2] has 8 free parameters, of which 5 parameters de- 
scribe the global cluster properties (n e0 , T e0 , r s , f3, n) 
and 3 additional parameters (r coo i, a, 7) are used to 
model the central region of the cool core clusters. 

3.3.1. Polytropic Clusters MS 1137.5+6625 and 
CL J 1226.9+3332 

For polytropic clusters, which do not have a cool core 
component, 5 parameters are sufficient to describe the 
distribution of density, temperature. We fixed /3 to 2 for 
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Fig. 3. — X-ray surface brightness and temperature profiles of MS 1137.5+6625 in the radial range 0-90 "and CL J1226. 9+3332 in the 
radial range - 120 ". The red line in both profiles shows the best fit model to the data; the green line in surface brightness profiles shows 
the background level. The overall x 2 °f the fit (Table [3j is the sum of the x 2 values of the surface brightness profile and the temperature 
profile. 

TABLE 3 
Best-fit parameters of the model 



Cluster 


™e0 
(10~ 2 cm- 3 ) 










T 
(keV) 






7 


X 2 (d.o.f.) 


P value 


MS 1137.5+6625 


^ -0.19 


23.84+l° i f 


4.83 + 


30 
27 


2.0 


7 84 + 1 - 08 
,,B -1 .26 








14.4 (26) 


96.7 % 


CL .11226.9+3332 


3 99+ ' 23 


22.65+_|'°f 


4.54+0 


54 
31 


2.0 


1 s qo + 2.28 








14.6 (29) 


98.8% 


Abell 2204 


4 42+ ' 37 
4 -* 2 -0.24 


ai.73±| : *S 


6.44+J 


02 
51 


1 39 + 04 
1 -0.06 


14.28l«; 7 = 


19 42 + - 60 
1J.42_ 73 


16+ 001 
ult> -0.01 


2.0 


115.5 (145) 


96.6% 


Abcll 1835 


„ -»+0.29 
-0.07 


40 - 32 ±6.48 


3.98 + JJ 


73 
41 


1 94+ ' 15 
l u -0.22 


i8. 2 8+; : « 


22 65 + ° 28 
22t>3 -1.17 


n ls +0.02 
018 -0.01 


2.0 


99.3 (93) 


30.8% 



these clusters (|Navarro "eTal][l996l) . since the polytropic 
index (n) and j3 cannot both be determined from X-ray 
data available. We report the x 2 of the best-fit model for 
MS 1137.5+6625 and CL J1226.9+3332 in Tabled We 
also calculate the gas mass by taking the volume integral 
of Equation IT31 and the total mass using Equation [2] and 
report the results in Table HJ 

3.3.2. Cool Core Clusters Abell 2204 and Abell 1835 

The clusters Abell 2204 and Abell 1835 have a clear 
cool core component (see Figure E} , which requires the 
use of the cooling equations described in ^2.41 The best 
fit model parameters are listed in Table [3] We also calcu- 



late the gas mass by taking the volume integral of Equa- 
tion [19] and the total mass using Equation [20l and report 
the results in Table @J 

4. COMPARISON WITH PREVIOUS WORK 

In Table [5] we present the comparison of mass mea- 
surements of Abell 1835 produced f rom the polytropic 
model with the masses reported in Mroczkowski ct al. 
(2009) at r 2 5oo and r 50Q . For this purpose we calcu- 
late gas mass and total mass at t he sa me radii r25oo 
and rsoo as in iMro czkows ki et all ()2009l ). and use the 
same Gaussian uncertainty on ta in order to have a fair 
comparison on masses. The gas mass measurements pro- 
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Fig. 4. — X-ray surface brightness and temperature profiles of cool core clusters Abell 2204 in the radial range - 540 "and Abell 1835 
in the radial range - 360 ". The red line in both profiles shows the best fit model to the data; the green line shows the background level. 
The overall \ 2 of the fit (Table [3]} is the sum of the x 2 values of the surface brightness profile and the temperature profile. 



TABLE 4 
Gas and total masses 



Cluster r 2 500 M ga s{r2&Oo) M*ot(r_50o) ^500 M gas (rsoa) M to t(r 5 oo) 
(arcsec) (10 13 M Q ) (10 14 M Q ) (arcsec) (1O 13 M ) (10 14 M Q ) 



MS 1137.5+6625 




i in+ - 13 
i - iu -o.n 


1 20+ - 36 
± ' zu -0.28 


98.5t!:§ 


3-08tri 3 3 


2 -ifi+ - 82 
z - oo -0.59 


CL J1226. 9+3332 


48.8±*:5 


o Qi+0.39 


2 16+ 069 




°- zy -0.63 


4 25+ 1 ' 22 
*- zo -1.14 


Abell 2204 


225.7±£1 


o QQ + 0.09 

j-y»_o.o9 


o o 7 +0.19 
— 0.18 


479.81^ 


10 S5+ 26 


fi 47+0-47 
0-*< _o.44 


Abell 1835 


150.6^1'^ 


^■ al -0.17 


o 79 +0.26 
J - ' -0.30 


309.7113 8 ! 


12.08t°; 3 « 


6 47+ 64 

"•^'-0.79 



duced by the polytropic model are consistent with the 
iMroczkowski et all (|2009f ) results at the la level (see Ta- 
ble ©. 

The Chandra Calibration Database (CALDB) has re- 
cently been revised to correct the effective area, result- 
ing in lower X-ray temperatures, especially for massive 
clusters 1 . The peak X- ray temperatures reported by 
IMroczkowski et all (|2009f ) using CALDB 3.4 are ~ 2 keV 
greater than the temperatures derived in this paper using 

1 http: / / cxc.harvard.edu/caldb/downloads/Release_notes / 
CALVB.v4AA.html 



the recent calibration (CALDB 4.1.1). From Equations[2] 
and[20l we estimate that this temper ature change would 
reduce the total masses reported by IMroczkowski et al.l 
(2009) by 17 % (see Table B). When the X-ray temper- 
ature calibration issue is accounted for, the total mass 
values in Table [5] are in agreement within the stated la 
uncertainties. 

5. DISCUSSION AND CONCLUSIONS 

We introduce a new model to describe the physical 
properties of the hot intra-cluster medium and present an 
application of the model to Chandra X-ray observations 
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TABLE 5 

Mass Comparison of Abell 1835 with IMroczkowski et al.I (120091) 





>"2500 

arcscc 


M2500, gas 

10 13 M ann 


M250O, tot 
10 14 M aun 


^500 

arcscc 


10 13 M sun 


M500,tot 

10 14 M aun 


Abell 1835 














Mroczkowski et al. (2009) 


I69.0±!;jj 


- -- + 0.25 
°- ' ' -0.35 


5.30±° f 2 [Note(a)] 


363.0 til " 


13.94t°;^ 


W.estl'lt [Note(b)] 


Polytropic Model (this work) 


169. Oj;™ 


5 79+ 035 

J - ' a -0.46 


4 13+O.31 

*- io — 0.32 


363.0 ±il;g 


14 31+ 071 

"■ ol -0 .82 


7 o 7 +0.82 
'•°'-0.83 



(a) When X-ray temperature is recalibrated using CALDB 4.1.1, mean value of A/2500. tot decreases to ~ 4.40 X 10 M s . 

(b) When X-ray temperature is recalibrated using CALDB 4.1.1, mean value of M500, tot decreases to ~ 8.86 X 10 14 M ev 



of MS 1137.5+6625, CL J1226.9+3332, Abell 1835 and 
Abell 2204. The model is based on a polytropic equation 
of state for the gas in hydrostatic equilibrium with the 
cluster gravitational potential. Using a function for the 
cluster total mass density that has the asymptotic slope 
as a free parameter, we obtain analytic expressions for 
the gas density, temperature and pressure. We also in- 
clude a core taper function that accounts for the cooling 
of the gas in the cluster center. 

This model has a number of features that make it suit- 
able for the analysis of X-ray and SZE observations of 
galaxy clusters. The model is analytic, and has a limited 
number of parameters which describe the global proper- 
ties of the cluster. For clusters which do not have a cool 
core, 5 parameters are sufficient to describe the distribu- 
tion of density, temperature, pressure and total matter 
density. The gas density and temperature are linked by 
the polytropic equation of state, and the total matter 
density is related to the plasma properties by the hydro- 
static equation. Therefore there is just one scale radius 
(r s ) that appears in the radial distribution of all thermo- 
dynamic quantities. The other parameters that describe 
the global physical properties of the cluster are the cen- 
tral density (n e o) and temperature (To) of the gas, the 
polytropic index n, and the asymptotic slope of the to- 
tal mass density (j3 + 1). For cool core clusters, three 



additional parameters allow an accurate description of 
the cooling of the gas in the core, and the accompanying 
increase in the density (§ 12.4)) . 

In addition to the analysis of spatially-resolved spec- 
troscopic and imaging X-ray data (see the model is 
applicable to SZE observations, which require a model 
for the plasma pressure. A number of models suit- 
able for SZE observation s are a vaila ble in the literature, 
for ex ample iNagai et al.I (|2007f ) and Mroczkowsk i et aLl 
(2009). Our model has the advantage of the simul- 
taneous applicability to both X-ray and SZE observa- 
tions, and it is therefore suitable for a number of cos- 
mological applica tions including the meas urement of the 
Hubble constant (|Bonamente et aI1 l2006'). the measure- 
ment of scaling relations betw een X-ray and SZE observ- 
ables (|Bonamente et al"1l2008l ). the measurement of clus- 
ter masses independent of cosmology from joint X-ray 
and SZE data (Hasler et al. 2010) and the measure- 
ment of the e ffect of He sedime ntation on X-ray mea- 
sured masses (jBulbul et al.lf2010h . 
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